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Abstract 

We consider a fourth order partial differential equation in n di- 
mensional space introduced by Abreu in the context of Kahler metrics 
on toric orbifolds. Similarity solutions depending only on the radial 
coordinate in R" are determined in terms of a second order ordinary 
differential equation. A local asymptotic analysis of solutions in the 
neighbourhood of singular points is carried out. The integrability (or 
otherwise) of Abreu's equation is discussed. 

1 Introduction 

In recent work Abreu has considered toric Kahler metrics on toric varieties 
or toric orbifolds ||5| of dimension 2n. Following a construction due to 
Guillemin each such variety or orbifold is completely determined by its 
moment polytope A in M", and the scalar curvature S of the Kahler metric 
is given by a formula 

S = -d.dug^^ (1.1) 

where the matrix elements g^^ are functions of x = (xi, . . . , x„)^ G M", and 
dj denotes the partial derivative d/dxj (the summation convention is as- 
sumed). More precisely, in terms of a potential function ^(x) = g{xi^ . . . , Xn) 
on M", the metric on the interior of the polytope A is 

ds^ = gjk dxjdxk, 

where gjk are the elements of the Hessian matrix G, i.e. 

G = (gjk), gjk = djdkg, 

and in ( |1.1D the curvature is determined in terms of the inverse matrix 
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In it is shown that the condition for extremal toric metrics is that the 
curvature is an affine function of x, in other words 



djS = constant, j = 1, . . . ,n. (1-2) 
Abreu constructs such metrics from potentials of the form 

1 

5 = -5^^Kx)log£Kx), (1.3) 
(=1 

where the ■^z(x) are affine functions which determine the facets of the poly- 
tope by the equations ■^/(x) = 0, / = l,...,d. For the simplest case of 
Guillemin's construction the polytope is a cuboid in M", so that d = 2n, 
and the toric variety is just (S'^)^^, i.e. n copies of the 2-sphere obtained 
by attaching an n-torus T"" to each point in the interior of the cuboid. 

The purpose of this note is to construct other types of solution to the 
equation (^]^) for the case of constant scalar curvature, namely 

djdkQ-^^ = —K, constant. (1-4) 



Henceforth we shall follow and refer to ( |1.4| ) as Abreu's equation, and 
show that for n > 1 it admits 0{n) invariant solutions of a different form 
compared with (pT^). 

Rather than being geometric, our motivation for considering the par- 
tial differential equation (PDE) (|1.4D comes from the theory of integrable 
systems. Taking derivatives of ( |l.3D we have 

where the constants cij are the coefficients of the affine functions ^i(x). 
The expression ( |l.5D is reminiscent of the potential for vanishing rational 
solutions |l^ of the integrable Kadomtsev-Petviashvili (KP) equation, 

dx \ dts dx dx^ J dt^ 

which takes the form 



w 



El dw , 



where t denotes the dependence of the poles on the times ^2)^3 (or an infinite 
sequence of such times in the full KP hierarchy). The dynamics of the d poles 
xiii) is governed by the integrable Calogero-Moser system for d particles. 
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The potential w in (1.7) is a central object in the Sato formulation of KP 



theory [|lq ] , since the pseudo-differential Lax operator L is constructed with 
a dressing operator W such that 

L = Wd^W-^, W = l+wd-^ + .... (1.8) 



The similarity between (|l.5| ) and the rational KP potential (|l.7D raises 
the question of whether Abreu's equation ( |1.4| ) might be integrable in some 
sense. Following the Ablowitz-Ramani-Segur conjecture ||2| that all reduc- 
tions of integrable PDEs should have the Painleve property, our natural 
instinct is to seek similarity reductions of ( |1.4| ) and examine the structure 
of their singularities in the complex plane. To avoid begging the question, 
by an integrable PDE we mean one having a Lax pair (ensuring solvability 
by the inverse scattering transform [^) and/or infinitely many symmetries 



[13|. 



2 Similarity solutions 



For each n, Abreu's equation (1.4) has 0{n) invariant solutions with the 
potential g being a function of the radial distance only, that is 

g = g[r), r=|x|. 

In that case the Hessian matrix and its inverse are given by 

r ~^ r (r) ' / ~'~ (r^/' r/) ' ^ 

where 

f = 9', 

and the prime denotes d/dr. 



Substituting the form ( |2.1| ) of the inverse Hessian matrix into (1.4) yields 
the following third order ordinary differential equation (ODE) for /: 

(" + rj^){^+rB' + („ + m)=-n: A = ^, B = ^ - ±^ 

(2.2) 

After an integration, this yields 

— + r5' + (n + l)S = Ar^"--, (2.3) 
r n 



for constant A. As a second order ODE for / = /(r), the equation (|2^ 
takes the explicit form 

V n ^ t I r 
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For the purposes of asymptotic analysis, (2.4) may be conveniently rewritten 



'^■^og(^) = (^r-Xr'-Ar. (2.5) 



as 

dr \ ^n-i J 

(By rescaling / we can always set k = 1, but we choose to leave k arbitrary.) 
In the case n = 1 it is straightforward to integrate ( |2.4D to obtain the general 
solution 

/ = - log ^ + a, 

p \p + X - Kr J 

so that the potential is 

g = />~^(r — r_) log(r — r_)+/j~"'^(r_(_ — r) log(r+— r)+ar+/5, r± = {X±p)/k, 

for arbitrary constants p > 0, a, (3. Up to shifting by an affine function, 
the solution for n = 1 is of the form ( |l.3D , with g{r) defined on a single 
interval (r_,r-|_) in M, which leads to a toric metric on 5^. Clearly for any 
n there is the trivial solution / = constant, for which the inverse Hessian 
in ( p.l| ) becomes infinite. For n > 1 we are unable to integrate ( |2.4| ) 
explicitly, and must resort to asymptotic analysis around singular points, 
before considering the solutions of the initial value problem. 
Asymptotics at r = 0: There are several different types of behaviour near 
the origin. 

• For any A, (^]^) admits the expansion 

/ ~ /o + ar" - 7^r"+i + 0{r^^\ (2.6) 
(n + 1) 

with arbitary constants /o 7^ and a. 

• For A 7^ there is an alternative expansion 

/ ~ ^0 + u ^ n ^""' + ^(''""'') (2.7) 
X[n — \) 

with /o 7^ 0, for n 7^ 2. When n = 2 we have instead 

/~/0 + A-V--i-r2 + O(r3). 
4AVo 



For A 7^ and n 7^ 2 there is another local expansion in (2.4) that is 
regular at r = 0, with leading order behaviour 



X[n — 1) 



for n 7^ 2. 



4 



In the case A = 0, there is an exact solution that is singular at the 
origin, namely 

/ = r . 

K 



Asymptotics at infinity: The equation {2A) admits two types of asymp- 
totic expansion at infinity: 

• In the first, / is asymptotic to a nonzero constant, i.e. 

/~/oo + ^^^r- + ..., (2.9) 
n 

for /oo ^ 0. 

• In the second case, / tends to zero: 

/ r-\ r^oo. (2.10) 

Movable singularities at r = tq / 0: The ODE (|2.4|) is outside the 
Painleve class of second order equations whose general solution has no mov- 
able singularities other than poles, as described in [^], since its solutions 
admit algebraic branching at arbitrary points rg 7^ in the complex r plane. 
The leading order behaviour in the neighbourhood of such a branch point is 



/~co(r-ro)i/", (2.11) 



where the constant cq is arbitrary. Since ( |2.4p is second order, we expect 
that ( 2.11 ) should be the first term in an expansion in powers of (r — ro)^/^ 



providing a local representation of the general solution, since it contains the 
two arbitrary constants rg and cq. If it has only algebraic branching around 
movable singularities, an ODE can possess the weak Painleve property as 



defined in |18|, and there are many examples of integrable ordinary and 
partial differential equations which have this property (see for instance 

HI). However, the ODE ( |2.4D also admits movable logarithmic branch 
points in its solutions, with 

in the neighbourhood of r = tq. Logarithmic branching such as ( p.l2| ) is 
taken as a strong indicator of non-integrability in differential equations . 

We would expect that generically the solutions of the ODE should have 
infinitely many branch points in the complex plane. If a solution has a 
branch point on the real axis, then for even n it cannot be real-valued for all 
real r. In that case, for A / it would make sense to consider a real- valued 
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branch defined on the interval [0, r*), where r* > is the position of the 
first positive real branch point. Such a solution might have the behaviour 
(2.8) or tend to a nonzero constant at the origin, and as r — > r* from below 



would look like 

/~c(r*-r)i/" (2.13) 

or 

/~clog(r*-r) (2.14) 

for real constants c, c. In that case the Hessian matrix G would be entirely 
determined by the function / and its first derivative, giving a metric on the 
ball of radius r* in M", with a singularity at r = r* where /' blows up. 
Exactly how to extend this to a metric on a 2n-dimensional (symplectic) 
manifold is not clear to us. 

In fact starting from initial data specified near the origin at r = e > 0, 
the ODE (|2.4|) is regular. It is not possible to take initial data at the origin 



since r = is a singular point of (2^). Clearly the right hand side of (|2^ 



is also singular at / = 0. It is easy to see that given /'(e) 7^ 0, the solution 
cannot have stationary points for r > e, since integrating (|2.4| ) with initial 
data /(e) 7^ 0, /'(e) = gives only the constant solution /(r) = /(e). It is 
then trivial to prove the following: 

Lemma: Suppose that initial data is specified for the ODE ( ^.4) cit « point 
r = e > 0, with /(e), /'(e) both nonzero. Then the solution f{r) for r > e 
remains positive (negative) for /(e) > (f{e) < 0), and for /'(e) > or 
/'(e) < it is monotone (increasing or decreasing, respectively), as long as 
it exists. 

We wish to consider solutions of the ODE which are free of branch points 
on the whole positive real axis, or on a finite interval [0, r*) (with branching 
at r = r*). In that case the components of the Hessian matrix G determine 
a metric on M", or on the ball of radius r*. The eigenvalues of G in ( |2.1| ) 
are //r (repeated n — 1 times) and /'. Thus if we also require a Riemannian 
metric given by a positive definite Hessian, then we should consider only 
solutions of (|2.4| ) defined with both initial data positive, i.e. /(e) > 0, 
/'(e) > 0. The signs of / and /' remain constant and the Lemma also holds 
when the solution is analytically continued to r < e, provided it exists. 
However, a priori we have no guarantee that when the solution is continued 
back towards r = (solving the initial value problem in reverse) it will 
not reach a singularity (branch point) at some point r = e' with < e' < 
e. Hence we are led to consider a nonlinear connection problem for ( |2.4| ), 
to determine what asymptotic behaviours at r = are compatible with a 
solution defined on [0,r*) with specified branching at r = r*, or a solution 
on the whole positive real axis with specified asymptotics at infinity. 

For completeness we will also consider the cases when either /, /' or 
both are negative, and henceforth we assume k > 0, A 7^ 0. Taking the four 
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different combinations of initial data in turn, we obtain the following result: 



Theorem 1: Suppose non-zero initial data is specified for the ODE ( 2.4 ) 
at r = € > 0. There are four possibilities: 

(i) f{e) > 0, /'(e) > 0. The solution reaches a logarithmic branch point with 
asymptotic behaviour ( 2. 1^ ) (with c < 0) at some finite point r = r* > e. 



(a) /(e) < 0, /'(e) > 0. The solution reaches an algebraic branch point with 
asymptotic behaviour (i2.l3i) (with c < 0) at some finite point r = r* > e. 
(Hi) /(e) > 0, /'(e) < 0. Either the solution reaches an algebraic branch 
point fi2.1$ ) (with c > 0) at some finite point r = r* > e, or f is defined 
for all r > e and asymptotes to a non-negative constant with the behaviour 
^2^) or ( \2.1(^ ) as r ^ oo. 

(iv) /(e) < 0, /'(e) < 0. Either the solution reaches a logarithmic branch 



point (2.14) (with c > 0) at some finite point r = r* > e, or f is defined for 
all r > e and asymptotes to a negative constant with the behaviour ( jg.^ j as 
r oo. 

Proof, (i) The solution cannot be defined for all r > e, since the only 
possible asymptotics (p.9|) or (2.10) at infinity necessitate /' < for suffi- 



ciently large r (with our assumption of positive k), which contradicts the 
Lemma. An algebraic branch point ( p.l3| ) would mean / ^ as r — > r*, 
contradicting the fact that / is monotone increasing by the Lemma. Hence 



the logarithmic branching (2.14) at some finite r = r* is the only possibility, 
(ii) Again by the Lemma, / must remain negative and monotone increasing. 
This rules out the asymptotic behaviours ( |2.9| ) or ( p. 10 ) at infinity for which 



/' < 0, and also a logarithmic branch point which would require / — > — cx). 
(iii) and (iv) are proved similarly. 

Having considered the standard initial value problem for (^^) at r = e, 
we can now solve it in reverse, and consider continuing the solutions of types 
(i)-(iv) backwards towards the origin. (This is equivalent to solving the ODE 
for /i(r) = /(— r) with initial data at r = — e.) The solution to the nonlinear 
connection problem may be summarized thus: 



Theorem 2: Solving the initial value problem for the ODE (2.4) in the 
reverse direction, with the four combinations of initial data as in Theorem 
1, leads to the following possibilities for the solution f{r) with r < e: 



(i) Either the solution reaches an algebraic branch point atr = r^ with ( \2.1l\ ) 
(cq > 0) for < ro < e, or it continues back to the origin with one of the 
asymptotic behaviours ( \2.(\ ) (with fo > 0, a > 0), ( \2. T[ ) (with /o > 0, for 
A > only) or (for A < 0, n / 2 only). 

(ii) Either the solution reaches a logarithmic branch point at r = rQ with 
( 2.12^ ) for < rQ < e, or it continues back to the origin with one of the 



asymptotic behaviours ( |g.6[ j (with fo < 0, a > 0), or 1{ ) (with /o < 0, for 
A > only). 
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fiii) Either the solution reaches a logarithmic branch point at r = rQ with 
) for < ro < e, or it continues back to the origin with one of the 
asymptotic behaviours ( jg. 6| j (with fo > 0, a < 0), or 1{) (with /q < 0, for 
A < only). 

(iv) Either the solution reaches an algebraic branch point at r = rQ with 
( 2.1i[ ) (Co < Oj for < ro < e, or it continues back to the origin with one 
of the asymptotic behaviours ( ^-dj ) (with fo < 0, a < 0), ([2. % ) (with fo < 0, 
for X <0 only) or (for A > 0, n / 2 only). 

Proof. The proof is straightforward, by a simple enumeration of the various 
possibihties ( |2.6D , ( |2.7| ) or ( |2.8D at r = allowed by the Lemma, as well as 
the different branching behaviours that can occur. 



3 Conclusions 

For any dimension n > 1, new solutions of Abreu's equation ( |1.4| ) have 
been found by making a similarity reduction to 0{n) invariant solutions. 
The similarity solutions are determined in terms of a single function /(r) 
depending only on the radial coordinate r, which satisfies a second order 
ODE (^^) including an arbitrary parameter A. Due to a singularity in 
the ODE at the origin, the initial value problem is not defined at r = 0. 
Nevertheless we have shown that for some suitable subset of the possible 
initial data specified at r = e > the solution may have a continuation in 
both the forward (r > e) and backward (r < e) directions to define a solution 
either on the whole positive real axis, or on an interval [0, r*) with a branch 
point at r = r*. With the restriction to Riemannian metrics defined by a 
positive definite Hessian matrix G in ( |2.lD , only the latter case is relevant 
(corresponding to case (i) in Theorems 1 and 2 above). However, whether 
such solutions could have further geometric significance, by extension to a 
suitable metric on a symplectic manifold of dimension 2n, is uncertain. 

Performing numerical integrations of (^^) with n = 3 and k = A = 1, by 
different choices of initial data we have obtained solutions displaying some 
of the asymptotic behaviours included in cases (i)-(iv) of the Theorems. A 
more detailed understanding of such solutions would necessitate estimates 
on the initial data to ensure that branch points do not appear as r ^ 0. 
Determining conditions for the existence of points of inflection would be an 
essential step in such an understanding. It would also be desirable to have 
some bounds on r*, the position of the branch point, in terms of /(e) and 
/'(e), but this is beyond the scope of the present work. 

The various Painleve tests are good heuristic tools for analysing ordi- 
nary and partial differential equations, and have been used to isolate new 
integrable systems (see Chapter 7 of for a review) . These tests are based 
on the Painleve property, that the general solution of an equation should 
have only poles as movable singularities. However, as originally empha- 



8 



sized by the authors of the Painleve property is extremely sensitive to 
changes of variables. For example, transformations of hodograph type can 
change equations with movable algebraic branch points into equations that 
have only poles |^ ^. The ODE (^3) has not only algebraic branching but 
also movable logarithmic branch points, and so certainly fails the standard 
Painleve test or its weak extension [O]. This is a strong indication that 



Abreu's equation (1.4) is not integrable. 

It has been observed |]l^ that when an equation is written in potential 
form it can have a single (or finitely many) logarithmic terms in a local 
expansion around singular points, and still be integrable. However, the 
asymptotic expansion for ( |2.4D with ( 2.12| ) as the first term seems to require 



infinitely many logarithms. In the case of the KP equation (|1.6D, the poten- 



tial w in (L8) is given by the logarithmic derivative of the tau- function: 



^ 1 

w = -— logr. 

ox 

The tau-function for the rational solutions (11.71) takes the form 



T = fl{x-xi{t)), (3.15) 
1=1 

and this is a polynomial both in x and the times t. For the more general 
algebro- geometric solutions of KP [11|, the tau-function is a theta- function 
of an arbitrary Riemann surface. The rational and soliton solutions arise 
as degenerate limits of the theta-functions. It would be interesting to see 
if Abreu's equation would admit quasiperiodic generalizations of (|l.5D , by 
dropping the extremality condition (|l]^) on the the curvature, allowing S 
in ( |l.ip to be a more general (say, periodic) function. 

Other methods of testing for integrability might provide useful informa- 



tion about Abreu's equation. Since (1.4) is not of evolution type, a promising 



method would be the symmetry approach of Shabat et al [Q3| , which has re- 



cently been extended [14| in order to deal with non-evolutionary equations. 
It would also be interesting to apply the methods of jl^ to look for other 
sorts of group invariant solutions of Abreu's equation ( |1.4| ). 
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